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ABSTRACT 

According to the second Fermi LAT Catalog ( 2 FGL), about one third of the 7 -ray sources listed 
have no assigned counterparts at lower energies. Many statistical methods have been developed to 
find proper counterparts for these sources. We explore the sky area covered at low radio frequency by 
Westerbork in the Southern Hemisphere (WISH) survey to search for blazar-like associations among 
the unidentified 7 -ray sources listed in the 2 FGL (UGSs). Searching the WISH and NRAO VLA Sky 
Survey (NVSS) radio surveys within the positional uncertainty regions of the 2 FGL UGSs, we select 
as 7 -ray blazar candidates the radio sources characterized by flat radio spectra between 352 MHz and 
1400 MHz. We propose new 7 -ray blazar associations for eight UGSs and we also discuss their spectral 
properties at low radio frequencies. We compare the radio flux density distribution of the low radio 
frequency 7 -ray blazar candidates with that of 7 -ray blazars associated with other methods. We find 
significant differences between these distributions. Finally, we discuss the results of this association 
method and its possible applicability to other regions of the sky and future radio surveys. 

Subject headings: galaxies: active - galaxies: BL Lacertae objects - radiation mechanisms: non-thermal 


1. INTRODUCTION 

In the recent decades, 7 -ray astrophysics has made 
great advances enabled by improvements in instrumenta¬ 
tion in the high energy technologies. The most recent and 
most accurate 7 -ray source catalog is the Fermi Large 
Area Telescope ( LAT) Second Source Catalog ( 2 FGL, 
iNolan et al.ll2012h . based on the first 24 months of data 
from LAT. Thanks to its silicon strip pair production 
and modern analysis processes, LAT has drastically re¬ 
duced the positional error of the sources with respect to 
previous studies, like those performed by the Energetic 
Gamma-Ray Experiment Telescope (EG RET ) on board 
the C ompton Gamma-Ray Observatory (jHartman et all 

mu. 

As the Fermi -LAT 7 -ray observations of the sky 
continue, the task of finding clear counterparts to 
all the detected sources becomes increasingly challeng¬ 
ing, mostly due to the large positional uncertainty of 
the faintest sources. Since a strong connection be¬ 
tween 7 -r ay and radio emi s sion h a s been clearly demon¬ 
strated (Ghirlanda e t all 120101: iMahonv et al.l 120101: 
lAckermann et al.l l 2011 at h it is natural to exploit radio 
surveys for finding new associations. Several associa¬ 
tion methods have been proposed to match the 7 -ray 
sourc es detected with s ource catalogs at low er frequen¬ 
cies (Pa ggi et al.l 2013; Masetti et al. I [20 1 3h . to give a 
proper counterpart to unidentified sources, where possi¬ 
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ble. Of the 1873 sources in the 2 FGL catalog, 575 have 
unknown physical nature. Regarding the Fermi sources 
of known types, blazars are the largest known population 
representing more than 80%. It is, therefore, fair to pre¬ 
sume that a significant fraction of the unidentified 7 -ray 
sources (UGSs) can be unrecognized blazar-like sources. 
Finding them is the aim of our work. 

Blazars are compact radio sources with flat (i.e., spec¬ 
tral index a < 0.5, where ~ u~ a ) radio spectra that 
steepen towards the infrared-optical bands; their overall 
spectral energy distribution shows two broadband com¬ 
ponents: the low energy one peaking in the infrared (IR) 
to X-ray frequency range, while the high energy one ex¬ 
tends up to the MeV-to-TeV range. Their emission fea¬ 
tures high and variable polarization, apparent superlumi¬ 
nal motions of radio jet components and high apparent 
luminosity, coupled with rapid flux variability over the 
whole electromagnetic spectrum. Blazars come into two 
classes: flat-spectrum radio quasars and BL Lac objects, 
which we label here as BZQs and BZBs resp ectively, 
following Roma-BZCAT (jMassaro et al.l l2QQ9h nomen¬ 
clature. Blazar emission interpreted as originating in 
radio loud active galaxies, with one of the relativistic 
jets emerging from the super-massive black hole in the 
galaxy nucleus pointi ng almost directly to the observer 
(jGiommi et al. (l 2012 f ). 

This paper is the latest of a seri es that focus on 
the nature of UGSs: iD’Abrusco et al. ( 201 3) (hereafter 
Paper I), which inves tigates sou rces as UGSs counter¬ 
parts as descr i bed in IMassaro et al.l (j2013al Paper II); 
IMassaro et al.l ( 2013 b. Paper III) adds another feature 
for the search of blazar-like sources f ocusing o n th e low- 
frequency radio feature of blazars; iPaggi et"ahl (j2013l . 
Paper IV) i nvolves t he X-ray emission as a distinctive 
feature and IMassaro et al ] d 2 om Paper V) propose a 
renewed IR approach, based on a 2-dimensional kernel 
density estimation (KDE) technique, all for the same 
purpose. 
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In thi s work, we apply the method proposed in 
[Paper III to search for blazar-like candidates within the 
7 -ray positional uncertainty regions of the UGSs listed in 
the 2 FGL and to select counterpart among them. This 
method is based on a multifrequency procedure, relying 
primarily on the flatness of the radio spectrum. The in¬ 
novative point in this method is that it is based on low- 
frequency (i.e., < 1 GHz) measurements, which is un¬ 
precedented in the investigation the e missio n of blazars 
and their association with UGSs. In [Paper III we com¬ 
bined the observations at 325 MHz from the Westerbork 
Northern Sky Survey (WENSS) at 325 MHz with those 
of the NRAO Very Large Array Sky Survey (NVSS) at 
1.4 GHz. In this work, we apply the same approach to 
the region covered by the counterpart of WENSS in the 
southern sky: the Westerbork in the Southern Hemi¬ 
sphere ("WISH. [De Brenck et al.ll2002[ ). Thanks to the 
combined results of this new search and of lPaoer III we 
build a reasonably sized population of new low-frequency 
selected blazar candidates. We study the spectral proper¬ 
ties of individual candidates and explore the flux density 
distribution of the whole sample in comparison to the 
y-ray blazars already listed in the 2nd F ermi LAT Cat¬ 
alog o f Active Galactic Nuclei f 2 LAC. fAckermann et al.l 
l2011b[) . We further refine our search by looking for IR, 
optical, UV and X-ray counterparts within the catalogs 
available. 

The paper is organized as follows: in Sect. 2 we de¬ 
scribe our method; in Sect. 3 we report and characterize 
our list of candidates; in Sect. 4 we discuss the spectral 
properties of all the sources found with the low-frequency 
method and make some projections regarding expecta¬ 
tions from new low-frequency radio telescopes and sur¬ 
veys. 

2. METHOD 

2 . 1 . Blazar characteristics in the low-frequency radio 

band 

The characteristics of blazars (in particular 7 -ray emit¬ 
ting blazars) in this ra nge of the electromagnetic spec¬ 
trum are discussed in iPaoer ILO : we briefly summarize 
those results here for the sake of clarity. 

Based on the cross correlation between the 2 FGL cat¬ 
alog, the Roma-BZCAT v4.1 (the most comprehensive 
blazar catalog up to now), and the WENSS, we de¬ 
fined one sample labeled as Low radio frequency Blazar 
(LB) and a subsample labeled Low radio frequency 7 -ray 
Blazar (LGB). Since the Roma-BZCAT catalog is based 
on the NVSS survey, we computed the low-frequency ra¬ 
dio spectral index values for these samples as: 

Q; - =log (^) /log (j_) (1) 

where v is the low radio frequency measured in MHz 
(i.e., 325 MHz for WENSS and 352 MHz for WISH), 
$1400 is the flux density measured at 1400 MHz in the 
NVSS survey and S v is the flux density measured at v 
frequency. Both flux densities are measured in mJy. 

The indices of about 80% of the sources from both the 
LB and LGB sam ples are smaller than 0 .5 and 99% have 
indices below 1.0 (Ma ssaro et a l. 2013b|). The flatness of 
the radio spectrum can be seen as a consequence of the 
dominance, also at low radio frequencies, of the inner jet 


departing from the core relative to the emission radiated 
by the larger structures of the blazar. 

In particular for the 7 -ray blazar sample, we consider 
as A class candidates the radio sources characterized 
by —1.0 < < 0.55 and B class as those with 

0.55 < < 0.65. These two classes have been de¬ 

fined to represent respectively the 80% ( A class ) and 
90% (B class ) of 7 -ray blazar radio spectral indices. In 
fact, the flatness of radio spectrum above 1 GHz is indeed 
a well known feature and was alread y used to discern 7 - 
ray blazar associations in the past (fHealev et al.l 12008 ). 
The low-frequency radio observations allow us to confirm 
this flatness down to 325 MHz. 

2.2. The WISH survey 

The WISH survey is the natural extension of the 
WENSS survey to 1.60 sr of the southern sky. at 
352 MHz. Both surveys were performed in the 1990s by 
the Westerbork Synthesis Radio Telescope (WSRT), us¬ 
ing the bandwidth synthesis mosaicing technique to com¬ 
bine 8 different bands of 5 MHz. WISH covers the area 
between — 9° < Dec < —26° to a limiting flux density 
of ^ 18 mJy (5cr), the same limiting flux density as the 
WENSS. Due to the low elevation of the observations, the 
survey has a lower resolution in declination (54" x csc(J)) 
than in right ascension (54"), resulting poorer source lo¬ 
calization than for WENSS by a factor ~ 2 x on average. 
Aside from this consideration, the WISH shares with the 
WENSS the same features upon which this method was 
calibrated, except for a negligible difference in frequency 
(An = 27 MHz). For this reason, we are confident in ap¬ 
plying the same association procedure that we used for 
WENSS. 

2.3. Association procedure 

For each UGS in the WISH footprint, we search for 
low-frequency sources in a circular region of radius equal 
to the semi-major axis of the 95% confidence level po¬ 
sitional uncertainty ellipse of the UGS itself. For each 
WISH source found, we look for a corresponding NVSS 
source in a circular region of radius equal to 8.5" and we 
then calculate the radio spectral index 0 ^ 352 ° • For this 
study, we only consider WISH sources classified as single 
component sources and featuring a NVSS association. 

In addition, we made local maps of the search regions 
for each UGS, overlaying on the WISH background the 
brightness contours of the NVSS map, the Fermi posi¬ 
tional uncertainty ellipse, and any possible blazar-like 
WISE detection up to 3.3" from the position of the 
matching NVSS source (see Fig. [Tj) . This gave us, in 
addition to a qualitative comparison of the relative posi¬ 
tions of various possible candidates and the UGS, a clue 
about their potential non-blazar nature if complex struc¬ 
tures are visible. Also, the angular separation between a 
candidate and the center of the Fermi ellipse was taken 
in account in cases of multiple matches, with the nearest 
source preferred. 

For associations between WISH and NVSS and be¬ 
tween NVSS and WISE, the search radii (of 8.5", and 
3.3", respectively) w ere selecte d based on a statistical 
study, as described in [Paper m. 

2.4. Multi-frequency data 
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Finally, we looked for additional multifrequency 
da ta for our new blazar candidates and those found 
in [Paper III In particular, we searched for matches 
with sources in the Australia T elescope 20 GHz sur- 
vey(AT20G, [Murphy et al.l 120101 1 for the WISH candi¬ 
dates, and in the_Green Bank 4.85 GHz northern sky 
survey (GB6, Gregor y et al.l[T996h for the WENSS can¬ 
didates. The AT20G survey, performed between 2005 
and 2006 with the Australia Telescope Compact Array 
(ATCA), contains sources with S 20.0 > 40 mJy; the GB6, 
performed between 1986 and 1987 with the NRAO seven- 
beam receiver, takes into account sources with £ 4.85 > 18 
mJy. For these surveys, source association is automat¬ 
ically provided by NASA/IPAC Extr agala ctic Database 
(NED). We also refer to iPaggi et al.l ( 2013 ) for the anal¬ 
ysis of X-ray emission observed by the Swift X-ray Tele¬ 
scope. 

3. RESULTS 

3.1. UGSs in the WISH footprint 

Of the 575 UGSs in the 2 FGL catalog, we discard all 
that feature a ‘c’ analysis flag to rule out potentially con¬ 
fused or spurious detections due to the strong interstellar 
emission ( INolan et al.ll 2012 f ). The resulting list of UGSs 
on the WISH footprint has 27 7 -ray sources. Of these, 
17 have at least one WISH and NVSS match (Table [l]). 

The total number of candidates is 31. Seven of these 
candidates have a radio spectral index in the range 
<0.55 and are considered class A candidates while 
six are in the range 0.55 < 0 ^ 25 ° < 0-65 and are consid¬ 
ered class B candidates. 

Sixteen out of the 31 candidates in the WISH foot¬ 
print feature a WISE detection in at least two IR bands; 
three of them are A class candidates and two are B class. 
However, we do not consider WISE detection to be nec¬ 
essary for selection. IR data (see Table [2]) are tak en from 
the AllWise November 2013 release [Cutri et al ] dnn). 
Flat spectrum radio quasars at high redshift or high syn¬ 
chrotron peaked BZBs could be too faint in the IR to 
be detected in one or more WISE bands; moreover, the 
relative astrometry of the WISH and NV SS might be 
poorer in this region (Massaro et al. 2013c) in respect to 
WENSS and NVSS. 

All considered, we propose eight UGSs blazar-like asso¬ 
ciations (i.e., five A and two B class candidates and one 
special case for WBN 2255.7—1833). In particular, we 
increase the number of class A and class B associations 
combined with [Paper Ill to 19 and 8 , respectively. 

3.2. Multiwavelength observations and radio flatness of 
low-frequency selected blazar candidates 

In this work and in [Paper Il3 . we evaluate the flatness 
of the radio spectrum according to Eqn. [l] on the basis 
of just two non-simultaneous flux density measurements, 
at 1.4 GHz (NVSS) and at low-frequency (352 MHz for 
the WISH or 325MHz for the WENSS). For this reason, 
it is important to find other radio flux density data to 
better investigate the spectral and variability properties 
of the candidates. We searched the literature for other 
radio surveys performed in the WENSS and WISH foot- 
prin ts; we consider here all the candidates found both in 
[Paper III and in the present work. 

In the WENSS footprint, we found that six A class 


and two B class candidates are listed in the GB 6 sur¬ 
vey at 5 GHz. We report them in Table [3] along with 
the radio spectral index, obtained with a weighted linear 
regression on the measurements at the three frequencies. 
For all these sources, the 5 GHz flux density is in good 
agreement with the extrapolation of the low-frequency 
spectrum, confirming that the radio spectrum is flat and 
flux density variability is not dramatic. A sample spec¬ 
trum is shown in Fig. [2] (left panel). 

In the WISH footprint, we find that one B class (WNB 
1251.8-2148, alias NVSS 125429-220419) is also listed 
in the AT 20 G. Data and radio spectral index regression 
are reported in Table [4] (see also Fig. [2j right panel). In 
this case, the extrapolation of the low-frequency power 
law clearly fails to match the high frequency data, sug¬ 
gesting prominent variability or a rather complex spec¬ 
tral shape, with a strongly inverted component above a 
few GHz; either way, the source behavior is consistent 
with the BZB class. 

4. DISCUSSION AND CONCLUSIONS 
4.1. Comparison with Paper III 

As we already noted, the WISH and WENSS surveys 
are very similar but, a posteriori, we can draw some use¬ 
ful conclusions regarding the impacts their different fea¬ 
tures had on our study. The poorer resolution in dec¬ 
lination of WISH (~ 2 x on average) mainly affects the 
spatial association and not the parameters on which the 
counterpart is selected. Since spatial association is cali¬ 
brated with a higher resolution survey like WENSS, ap¬ 
plying it to the WISH region and catalog results in a 
more conservative approach. The fact that we found al¬ 
most the same ratio (i.e. ^ 1 . 8 ) of candidates per UGS in 
both WENSS and WISH (58 WENSS candidates for 32 
UGSs and 31 WISH candidates for 17 UGSs) indicates 
that, on average, spatial association is not compromised. 

Similarly, the resolution discrepancy cannot be blamed 
for any impact on the quality of the associations pro¬ 
posed among the candidates, given the similar rate of 
7 -ray blazar -like associations in the two surveys. In fact, 
in [Paper IIJ we propose 21 (including one spe cial associ - 
ation that is neither class A nor class B\ see Pap er Ilj ) 
new 7 -ray blazars out of 65 UGSs in the WENSS region 
(i.e., 32%) and 8 out of 27 in WISH region (i.e., 30%)Q. 

4.2. Comparison with other methods 

Other methods have been suggested to recognize low- 
energy counterparts to UGSs, or at least to provide a 
statistica lly significant class i ficatio n of these sources. For 
instance, lAckermann et al.l ( 2012 ) have developed a sta¬ 
tistical approach to clas sify UGSs in the first catalog 
of Fermi sources ( 1 FGL, lAbdo et aD 112010) . Six of the 
2 FGL UGSs in the WISH footprint are associated with 
1FGL sources analyzed bv lAckermann et al.l ( 2 012 l h five 
of which are AGN-like and one of which is pulsar-like. 
Within the former sample, there are two sources for 
which we propose an association with a blazar counter¬ 
part on the basis of the low-frequency spectrum: 2 FGL 

7 The small discrepancy [with respect to what?] can be ex¬ 
plained, in addition to simple stati stics, on the grounds of the 
slightly different starting list, since in [Paper IIJ we considered only 
7 -ray sources without any analysis flag, while here we excluded only 
the sources with the c 7 -ray analysis flag in the 2FGL (INolan et al.l 
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J2017.5—1618 and 2 FGL J235 8.4-1811 . Intere stingly 
for the single UGS for which 1 Ackermann et alJ ( 201 2 ) 
propose a pulsar classification ( 2 FGL J1544.5—1126), our 
method finds a WENSS-NVSS match with quite steep 
spectral index (a = 0.74 =b 0.03), rejecting a blazar sce¬ 
nario. 

We further compared our results with the proposed 
associations found in the other papers of this series. In 
Figure [3] we show the comparison between the distribu¬ 
tion in the IR [3.4] — [4.6] — [12] fi m color-color plane 
provided by WISE for the 7 -ray emitting blazars and 
the sources selected in this work. The overall distribu¬ 
tion of the whole set of the simultaneous WISH-NVSS 
matches (black, red and green dots) is quite more scat¬ 
tered than the 7 -ray blazar population (orange dots). 
However, when we only consider the low-frequency se¬ 
lected blazar candidates, all of the most prominent out¬ 
liers are excluded and the remaining five sources (black 
and red dots) are in much better agreement with the IR 
colors of 7 -ray blazars. In two cases, the A-class candi¬ 
dates NVSS 120900-231335 and NVSS 222830-163643, 
the agreement is very good. For the latter source, associ¬ 
ated with 2 FGL J2228.6—1633, our method also provides 
the same counterpart candidate selected on the basis of 
the kernel density estimator techniq ue to IR col ors o f 
WISE cou nterparts applied to X -ray (jPaggi et al.l l2013h 
and radio (Massaro et al. 2013c) data. 

4.3. Radio flux density analysis 

The 7 -ray blazar candidates selected in this work and 
in llPap er H3 have by selection the same radio spectral 
properties of confirmed Roma-BZCAT blazars and of 
2 FGL blazar associations. It is natural to wonder why 
these sources have not been detected in the 2 FGL cata¬ 
log and eventually associated, e.g. in the second catalo g 
of AGNs detected by Fermi (jAckermann et al.l l 2011 bh . 
The sources could have been excluded from the 2LAC 
because they do not formally pass the threshold for being 
considered high confidence associations (a test that does 
not take into account the spectral index); typically, this 
could happen for low flux density radio sources, which 
have a larger spatial density. It is thus likely that, to¬ 
gether with a general similarity to the already known 
blazars, our candidates also have some peculiar charac¬ 
teristics. 

For this reason, we show in Fig. 0] the distribution 
of radio flux density at 1.4 GHz for all the blazars in 
the 2 LAC and for the candidates selected here and in 
iPaoer III The two distributions are clearly different, as 
confirmed by a KolmogorovSmirnov test which yields 
a probability of 2.7 x 10 -15 of being obtained from 
the same population. In particular, the 2LAC blazar 
flux density distribution is shifted to much larger val¬ 
ues. This strongly suggests that our method is very ef¬ 
ficient at selecting faint blazars. These sources are po¬ 
tentially of great interest: if they are of the BZB type, it 
could mean that they could be of the extreme and elu¬ 
sive class of ultra-high synchrotron peaked sources; one 
prominent example could be WNB 2225.8—1652, which 
has an inverted low-frequency radio spectral index of 
a = —0.19 4= 0.10 and is our proposed association for 
the UGS 2 FGL J2228.6—1633, characterized by a 7 -ray 
spectrum as hard as T = 2.07 ±0.16. On the other hand, 
some of our candidates could also be faint BZQs, and in 


this case their low flux densities could stem from their 
high redshifts (like WN 1500.0±4815 with an estimated 
redshift of 2.78, see Mass aro et al.ll2013bf h which would 
also make them of scientific value. No WISH sources 
have a measured redshift. 

4.4. Summary and outlook 

We have searched for counterparts of the the 27 UGSs 
in the 1.6 sr footprint of theyWISH survey, following the 
methods described in iPaper hi and based on the flat 
spectrum at low-frequency characteristic of blazars. We 
have found blazar-like associations for 8 UGSs, that to¬ 
gether with the 23 sources selected in the WENSS foot¬ 
print with the same method provides a blazar association 
for a sample of 30 new 7 -ray blazar associations. This 
sample extends the distribution of the radio flux den¬ 
sity of 7 -ray blazars to lower values, allowing us to study 
otherwise elusive AGNs. 

The application of our method thus shows promis¬ 
ing results in terms of numbers of counterparts pro¬ 
posed and their physical features. In particular, the 
possibility of using low-frequency radio data to find 
UGSs counterparts is even more important in the light 
of the imminent start of numerous studies in this fre¬ 
quency range, like the L Ow Frequency ARray (LOFAR, 
Ivan Haarlem et aLll2013lh the M urchison Widefield Ar¬ 
ray fMWA. iTingav et al l 120131) , the Long Wavelength 
Array (LWA, Elling son et al.l 120091 ). and eventually 
the S quare Kilometer Array (SKA, e.g. iDewdnev et al.l 
ESmJ). These facilities will allow our method to be ex¬ 
tended using an even deeper and simultaneous dataset, 
while dedicated targeting of the candidates with opti¬ 
cal spectroscopy and Very Long Baseline Interferometry 
observations will confirm the natures of the proposed as¬ 
sociations. 
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APPENDIX 

Notes on individual UGSs 

Below, we report a brief analysis of the results for each UGS listed in Table [lj 

• 2FGL J0340.7— 2421 Both candidates in this search region have an IR detection by WISE. Our method selects 
WN 0338.4—2436 as a blazar association due to its flatter ( class A) spectral index; this source also has a greater 
flux density. 

• 2FGL J0600.8—1949 There are three candidates: WN 0558.8—1950 provides the best association among the 
others since it is nearest to the 7 -ray position and it has blazar-like IR colors, a flatter spectrum, and larger flux 
density; however, the radio spectrum is steeper than 0.71 so we can not formally classify it as a blazar. 

• 2FGL J1059.9—2051 Both candidates in the search region are detected in the IR by WISE&nd have large flux 
density but have a quite steep radio spectral index. 

• 2FGL J1208.6—2257 WN 1206.4—2256 is the best association among the other candidates because of its flat 
(class A) spectral index and WISE blazar-like detection. 

• 2FGL J1254-2—2203 WNB 1251.8—2148 has an intermediate spectral index and is proposed as a class B blazar 
association. We note that a second nearby 1.4 GHz source (NVSS 125422-220413, see Fig. [lj right panel) is not 
present in the WISH survey, suggesting a flatter spectral index; this source could also be a blazar-like candidate 
and a partial contributor to the 7 -ray emission in this region. 

• 2FGL J1458.5—2121 The only candidate has the IR features of a blazar-like source but a too-steep spectral 
index; moreover, a much brighter source that could be the association for the UGS lies just outside the search 
region. 

• 2FGL J1544-5—H26 The only candidate has a too-steep radio spectrum and no proper WISE match, so it is 
n ot a l i kely blazar a s soci ation. Indeed, this UGS was already detected in the 1 FGL and the statistical method 
of I Ackermann et~ahl ( 20121 ) suggested that it should be classified as a pulsar. 

• 2FGL J1624-2—2124 Four candidates are present within the 7 -ray error radius, some with A or B class spectral 
index, and one with a WISE detection. The most promising potential associations are WNB 1620.7—2120 (flattest 
a) and WNB 1621.1—2119 (largest flux density and presence of IR emission); we prefer the first, consistent with 
the method. 

• 2FGL J1631.0—1050 None of the three candidates within the search region has a flat spectral index. 
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Fig. 1.— Map of the search regions around 2FGL J1208—2257 (left panel) and 2FGL J1254.2—2203 (right panel). Color scale background 
shows the WISE image, contours represent the 352 MHz emission from the WISH, and crosses indicate NVSS sources. The white ellipse is 
the 7 -ray 95% confidence region. In the right-hnad panel, our candidate is the leftmost of the twin NVSS sources. (Note that the righthand 
source is formally not in the WISH catalog.) 
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Fig. 2 .— Radio spectra of NVSS J030727+491510 (left panel) and NVSS J125429—220419 (right panel). The solid lines represent the 
linear regression spectral index, which is clearly a good fit to the data only in the first case. 

• 2FGL J1646.7—1 333 WN 1644.0—1323, even with a WISE detection and typical blazar-like X emission 
(|Paggi et alil2013 ), has a too-steep radio spectral index to be considered a blazar-like association. 

• 2FGL J1913.8—1237 According to our method, none of the candidates in this region can be considered a good 
blazar association. 

• 2FGL J2009.2—1505 The only candidate, even featuring blazar-like X emission (|Paggi et al.lj201 3). has a some¬ 
what too steep spectral index for being considered a blazar-like source, and it also lacks a WISE detection. 

• 2FGL J2017.5—1618 WNB 2014.9 —16 2 7 is a class A blazar association, even if it would not be selected as such 
on the basis of the IR emission. lAckermann et al ] (j2012h classify this UGS as a likely AGN. 

• 2FGL J2031.4~ 1842 Both candidates in the search region have a class B radio spectral index; the most likely 
blazar association is WNB 2027.8—1900, which is brighter and detected in the IR. 

• 2FGL J2124-0—1513 The only candidate has a too-steep spectral index for being selected as a blazar-like source, 
even though it has a WISE detection. 

• 2FGL J2 228.6— 1633 The inverted spectral index of WNB 2225.8—1652, its blazar-like WISE colors, and X-ray 
emission (jPaggi et al.ll2013f ) make this source a highly reliable blazar association for the UGS. 

• 2FGL J2358.4—1811 WNB 2355.7—1833 has a radio spectral index just above our B class threshold and a WISE 
detection. Although its IR colors are not very blazar-like, the error range on the spectral index overlaps with 
the B class range. We consider it a suitable, even if not forma l, blazar-like association, in agreement with the 
AGN classification proposed for this UGS bv lAckermann et al. ( 201 2). 
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Fig. 3. — The [3.4] — [4.6] — [12] jtmi color-color plot of all the WISE counterparts of the WISH-NVSS sources (green diamonds) and 
candidate blazars (black squares: class A; red squares: class B ) in comparison with the blazars that constitute the WISE 7 -ray strip (orange 
dots). 




Log F (mJy) Log f 


Fig. 4.— Normalized distributions of flux density at 1.4 GHz. Left panel: 7 -ray blazar candidates selected among the UGSs in the 
WENSS and WISH footprints; right panel: sources in the 2LAC (solid line: all 2LAC sources, dotted line: BZBs, dashed line: BZQs). 
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TABLE 1 

Possible blazar-like counterparts for UGSs in the WISH region. For each source, we list the 2FGL name (Col. 1), the list 
OF CANDIDATES PER UGS WITH WISH NAME (COL. 2) AND NVSS NAME (COL. 3), THE CORRESPONDING FLUX DENSITY AT 352 MHZ (COL. 
4) and 1400 MHz (Col. 5), the presence in the WISE catalog (Col. 6), the low-frequency radio spectral index (Col. 7) and 

THE RELATIVE CLASSIFICATION (COL. 8). 


2FGL name 

WISH name 

NVSS name 

S352 

(mJy) 

£1400 

(mJy) 

WISE 

detection 

M400 

a 352 

Index 

Class 

2FGL J0340.7—2421 

WNB 0338.0-2425 

NVSS 034011-241602 

72 ±4 

23.5 ±1.1 

V 

0.81 ±0.04 



WNB 0338.4-2436 

NVSS 034033-242712 

111 ±4 

82 ±3 

v 

0.22 ± 0.03 

A 

2FGL J0600.8—1949 

WNB 0557.9-1954 

NVSS 060003-195446 

44 ±4 

12.7 ±0.6 

V 

0.90 ±0.06 



WNB 0558.8-1950 

NVSS 060100-195049 

256 ±4 

96 ±3 

V 

0.71 ±0.02 



WNB 0559.4-1948 

NVSS 060138-194853 

124 ±4 

43.2 ±1.4 


0.76 ± 0.03 


2FGL J1059.9-2051 

WNB 1057.1-2037 

NVSS 105935-205311 

414 ±8 

138 ±5 

V 

0.80 ± 0.02 



WNB 1057.7-2040 

NVSS 110014-205621 

737 ±4 

232 ±7 

V 

0.84 ±0.02 


2FGL J1208.6—2257 

WNB 1205.6-2232 

NVSS 120816-224925 

151 ±3 

66 ±2 


0.60 ± 0.02 

B 


WNB 1206.4-2256 

NVSS 120900-231335 

34 ±3 

28.8 ± 1.0 

V 

0.12 ±0.06 

A 


WNB 1207.5-2234 

NVSS 121007-225106 

96 ±3 

35 ± 1.1 


0.73 ± 0.03 


2FGL J1254.2—2203 

WNB 1251.8-2148 

NVSS 125429-220419 

40 ±4 

18.4 ±0.7 


0.56 ± 0.06 

B 

2FGL J1458.5—2121 

WNB 1456.2-2112 

NVSS 145904-212357 

190 ±3 

70 ±2 

V 

0.73 ± 0.02 


2FGL J1544.5—1126 

WNB 1541.4-1115 

NVSS 154414-112443 

83 ±3 

30.0 ± 1.0 


0.74 ± 0.03 


2FGL J1624.2—2124 

WNB 1620.5-2058 

NVSS 162332-210457 

19 ±3 

9.9 ±0.6 


0.47 ±0.10 

A 


WNB 1620.7-2120 

NVSS 162345-212716 

30 ±3 

19.2 ±1.5 


0.32 ± 0.07 

A 


WNB 1621.1-2119 

NVSS 162403-212645 

175 ±3 

80 ±3 

V 

0.57 ±0.02 

B 


WNB 1623.0-2111 

NVSS 162600-211825 

105 ±3 

45.2 ± 1.4 


0.61 ±0.03 

B 

2FGL J1631.0-1050 

WNB 1628.0-1045 

NVSS 163049-105218 

106 ±5 

29.7 ± 1.0 

V 

0.92 ± 0.03 



WNB 1628.7-1037 

NVSS 163130-104322 

129 ±5 

45.6 ±1.8 


0.75 ± 0.03 



WNB 1628.8-1044 

NVSS 163139-105057 

591 ±5 

184 ±6 

V 

0.85 ± 0.02 


2FGL J1646.7—1333 

WNB 1644.0-1323 

NVSS 164651-132849 

447 ±4 

99 ±3 

V 

1.09 ± 0.02 


2FGL J1913.8-1237 

WNB 1910.5-1235 

NVSS 191320-122949 

44 ±4 

14.7 ±0.6 


0.79 ± 0.06 



WNB 1910.8-1246 

NVSS 191339-124120 

133 ±4 

39.7 ± 1.3 


0.88 ± 0.03 


2FGL J2009.2-1505 

WNB 2005.8-1513 

NVSS 200838-150500 

299 ±6 

105 ±4 


0.76 ± 0.02 


2FGL J2017.5—1618 

WNB 2014.9-1627 

NVSS 201745-161820 

110 ±3 

53.6 ± 1.7 


0.52 ± 0.02 

A 

2FGL J2031.4—1842 

WNB 2027.6-1850 

NVSS 203030-184033 

204 ±3 

92 ±3 


0.58 ± 0.02 

B 


WNB 2027.8-1900 

NVSS 203044-185033 

324 ±3 

134 ±4 

V 

0.64 ± 0.02 

B 

2FGL J2124.0—1513 

WNB 2121.8-1533 

NVSS 212438-152017 

147 ±3 

40.0 ± 1.3 

V 

0.94 ± 0.02 


2FGL J2228.6-1633 

WNB 2225.8-1652 

NVSS 222830-163643 

16 ±3 

20.9 ± 1.1 

V 

-0.19 ±0.10 

A 


WNB 2226.0-1641 

NVSS 222842-162619 

16 ±3 

7.2 ±0.5 


0.58 ±0.10 

A 

2FGL J2358.4—1811 

WNB 2355.7-1833 

NVSS 235820-181621 

25 ±3 

9.7 ±0.6 

V 

0.69 ± 0.09 



TABLE 2 

WISE [3.4]-[4.6] and [4.6]-[12] magnitudes. Data are expressed in terms of magnitudes (Vega system) and represent the 

DIFFERENCE BETWEEN TOTAL IN-BAND BRIGHTNESS MEASUREMENTS. THE [4.6]-[12] ERROR IS NOT PRESENT WHENEVER THE [12] BAND 

DETECTION IS AN UPPER LIMIT. 

NVSS name 

[3.4]-[4.6] mag 

[3.4]-[4.6] error mag 

[4.6]-[12] mag 

[4.6]-[12] error mag 

NVSS 034011-241602 

0.59 

0.23 

3.81 


NVSS 034033-242712 

0.23 

0.08 

2.73 


NVSS 060003-195446 

0.44 

0.08 

3.30 


NVSS 060100-195049 

0.05 

0.09 

2.83 


NVSS 105935-205311 

0.74 

0.16 

3.04 


NVSS 110014-205621 

0.39 

0.27 

4.19 


NVSS 120900-231335 

1.06 

0.04 

3.27 

0.08 

NVSS 145904-212357 

1.15 

0.22 

3.69 


NVSS 162403-212645 

0.26 

0.07 

1.56 


NVSS 163049-105218 

0.70 

0.44 

4.95 


NVSS 163139-105057 

0.00 

0.05 

1.14 


NVSS 164651-132849 

0.67 

0.04 

2.11 

0.19 

NVSS 203044-185033 

0.78 

0.11 

3.01 

0.38 

NVSS 212438-152017 

0.38 

0.44 

4.91 


NVSS 222830-163643 

0.75 

0.04 

2.37 

0.14 

NVSS 235820-181621 

0.18 

0.15 

3.42 
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TABLE 3 

Regression of radio spectral index with GB6 data of WENSS possible blazar-like counterparts. 


NVSS name 

Class 

log $325 

[mJy] 

log $1400 

[mJy] 

log $4850 

[mJy] 

a 

regr 

^1400 

a 325 

NVSS J030727+491510 

A 

1.90 

± 

0.05 

1.75 

± 

0.03 

1.53 

± 

0.15 

0.26 

± 

0.08 

0.24 ± 

0.04 

NVSS J033153+630814 

A 

1.79 

± 

0.05 

1.63 

± 

0.03 

1.40 

± 

0.16 

0.28 

± 

0.09 

0.26 ± 

0.04 

NVSS J035309+565431 

A 

1.98 

± 

0.05 

1.76 

± 

0.03 

1.64 

± 

0.11 

0.32 

± 

0.08 

0.34 ± 

0.04 

NVSS J072354+285930 

A 

1.90 

± 

0.06 

1.56 

± 

0.03 

1.49 

± 

0.16 

0.48 

± 

0.10 

0.53 ± 

0.04 

NVSS J150229+555204 

A 

1.70 

± 

0.06 

1.54 

± 

0.03 

1.36 

± 

0.17 

0.26 

± 

0.10 

0.25 ± 

0.05 

NVSS J150147+480335 

A 

1.43 

± 

0.12 

1.32 

± 

0.03 

1.41 

± 

0.15 

0.06 

± 

0.16 

0.18 ± 

0.08 

NVSS J210805+365526 

A 

1.83 

± 

0.05 

1.88 

± 

0.03 

1.83 

± 

0.10 

-0.04 ± 0.08 

-0.08 ± 0.04 

NVSS J060102+383828 

B 

3.261 

± 

0.002 

2.85 

± 

0.03 

2.51 

± 

0.09 

0.65 

± 

0.48 

0.65 ± 

0.02 

NVSS J101657+560112 

B 

2.00 

± 

0.04 

1.61 

± 

0.03 

1.79 

± 

0.10 

0.44 

± 

0.06 

0.62 ± 

0.03 


TABLE 4 

Regression of radio spectral index with AT20G data of WISH possible blazar-like counterparts. 


NVSS name 

Class 

log $352 [mJy] 

log $1400 [mJy] 

log $20000 [mJy] 

QTegr 

n4400 

a 352 

NVSS J125429—220419 

B 

1.60 ±0.10 

1.26 ± 0.04 

1.85 ±0.07 

-0.31 ±0.06 

0.56 ± 0.06 



